Theory of anomalous magnetic interference pattern in mesoscopic 
superconducting/normal/superconducting Josephson junctions 
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The magnetic interference pattern in mesoscopic SNS Josephson junctions is sensitive to the 
scattering in the normal part of the system. In this paper we investigate it, generalizing Ishii's 
formula for current-phase dependence to the case of normal scattering at NS boundaries in an SNS 
junction of finite width. The resulting flattening of the first diffraction peak is consistent with 
experimental data for S-2DEG-S mesoscopic junctions. 



I. INTRODUCTION 

The Fraunhofer diffraction pattern characterizing the 
dependence of the critical Josephson current on an 
applied external magnetic flux is a well-known fea- 
ture of tunneling or dirty superconducting-normal- 
superconducting (SNS) junctionsi. The latter is defined 
by the condition that quasiparticle scattering length by 
impurities in the normal layer, is much less than the 
characteristic dimensions {W, L) of this layer. For the 
case of clean SNS junctions, this pattern is more complex 
and demonstrates a distinct triangular central peakS. 
The pattern is <i>o-periodic in either case. More recent 
experimental and theoretical efforts have revealed, that 
in mesoscopic SNS junctions this pattern may be strongly 
altered by geometrical effects. For example, Heida et 
al- observed a 2$o-periodic diffraction pattern in narrow 
junctions, in contrast to the expected $o- The essential 
features of this behavior were explained theoretically^ 
using a semiclassical technique that accounts for the re- 
stricted geometry in terms of the classical trajectories 
associated with current-carrying Andreev states states 
in the normal layer. 

In this paper we generalize this approach to include the 
effects of normal reflection at the NS boundaries, which 
are treated in the Blonder-Tinkham-Klapwijk (BTK) ap- 
proximation^. The motivation for this work was provided 
by the JHAT (Jensen, Harada, Akazaki, and Takayanagi) 
experiment^. They observed critical current oscillations 
vs. flux in a wide ballistic S-2DEG-S junction, which ex- 
hibited apparently random small-scale (a fraction of $o) 
oscillations on top of a flattened flrst interference peak. 
As the flrst step towards understanding this data, here 
we extend previous work—'- (which has focused mainly on 
the case without normal reflection at the NS boundary) 
by including closed trajectories formed by reflection from 
the walls of the device as well as from the NS boundaries. 
Such trajectories may be important because, for them, 
the effective area penetrated by the magnetic field (and 
therefore the field-dependent phase of the corresponding 
contribution to the supercurrent) change smoothly with 
their variations. By contrast, the contribution due to 
chaotic trajectories is expected to exhibit a rapid change 



with variation of the trajectory and thus the field depen- 
dence effectively cancels on scales below $0i like in an 
SNS junction in the diffusive limit. We show, that such 
an approach captures some, but not all, essential features 
of the experiment. 

The general idea behind the quasiclassical approach is 
to represent the Josephson current through the SNS junc- 
tion as a sum of contributions from various quasiclassical 
trajectories linking the superconductors'*. Such a gener- 
alization of Ishii's formulck^ follows from the quasiclassical 
Eilenberger equations and gives a transparent descrip- 
tion suitable for complicated geometries. We will assume 
that the conditions of its apphcabihty (Af ^ *C M^, L, 
where Af is the Fermi wavelength and is the supercon- 
ducting coherence length ["Cooper pair size"]) are satis- 
fied. Our approach is to first calculate the density N{£^) 
of Andreev states in the normal area and then express the 
Josephson current through it. (Here ^ is the quasiparticle 
energy measured from the Fermi level.) For simplicity, 
we consider the case of a rectangular normal region with 
dimensions W x L (see Fig.^, with specular normal re- 
flections from the boundaries (given by the line segments 
CD and FE in Fig.^. 

The flux (f> — HWL produced by the magnetic fleld 
H = Hcz through the normal region can be described by 
the vector potential 



A = -Hy 



(1) 



The applied magnetic field is screened by the super- 
conducting regions, producing screening currents over a 
length given by A, the magnetic penetration depth. As is 
standard* , we shall neglect the variation of the supercon- 
ducting phase and magnetic field over this length scale; 
thus, our results apply at W,L ^ A. 

Before turning to the calculation of the current, let us 
briefiy remark on why closed trajectories are of princi- 
pal importance. The introduction of normal reflection 
on NS boundary takes away the simplicity of quasiclassi- 
cal trajectories with Andreev reflections only (when the 
hole retraces the electron's path). Even if assume (as we 
do) , that the Andreev reflection exactly reverts the group 
velocity of a quasiparticle, and that a normal reflection 
is specular, the quasiparticle trajectory will be generally 
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FIG. 1: Schematic picture of a closed quasiparticle trajectory 
carrying supercurrent between the superconducting leads (la- 
belled S). The dashed line is at i/ = 0. 



w and w along AOB and BQA, via 
w • $0 = / A -dl =-H 



A-dl 

AOB 

H{Saobef 



y{x)dx, 



w ■ $0 



WL/2), (2) 
A - dl ^ -H / y{x)dx, 

BQA J L 

-H{Sbqafe - WL/2). (3) 



We have expressed the appropriate area in Eq. ^ [and 
similarly in Eq. Q] as the difference of the area below 
AOB (which we denote by Saobef) a-nd the area below 
the dashed line. Our expression for the Josephson current 
carried by a particular trajectory will involve both the 
sum, V = w + and difference, = it; — w, of these 
phases: 

1/ • $0 = —H{Saobef — Sbqafe), (4) 

= -HSaOBQ = -$e:r, (5) 

17-$o = -H{Saobef + Sbqafe) + HWL, (6) 



open, with multiple reflection points at an NS boundary. 
Each SNS leg of such a trajectory coherently contributes 
to the latter's partial Josephson current Itr., with wildly 
varying phase factors. Therefore the total contribution 
of such trajectories will be suppressed. The exceptions 
are closed trajectories, where phase gains are systematic. 
Such trajectories contribute most to the quasiclassical 
density of stateaSi. If normal reflection amplitude from 
NS boundary is small, \TZ\ <^ 1, the leading contribution 
will be from the trajectories with the minimum possible 
number of reflections on the NS boundaries (one each, 
see Fig. ^ (class A trajectories) , which simplifies the sit- 
uation. On the other hand, in this case we can not dis- 
card the contributions from open (class B) trajectories, 
since only the first leg will significantly contribute to the 
current, and the above "phase cancelling" argument no 
longer applies. 



In the following, we assume that 1) Andreev reflection 
on NS boundary prevails and 2) the scattering in the nor- 
mal region is negligiblaiS. We will therefore concentrate 
on the two above classes of trajectories. Consider then 
the closed trajectory AOBQ in Fig. ^ where electrons 
propagate clockwise, and holes, related to the electrons 
by Andreev reflections from NS boundaries (at A and B 
in Fig.nj, propagate counterclockwise. (The results pre- 
sented here are valid for a trajectory with an arbitrary 
number of reflections from the side walls). The phase 
gain for an electron from the magnetic field along AOB 
is simply related to the flux enclosed by the trajectory 
and the dashed line y = 0. We introduce the phase gains 



where ^ex is the flux of the external field through the 
loop Saobq- Let us compute the phase gains (due to 
the field and from propagation) accumulated by quasi- 
particles travelling along AOB and AQB. Denoting the 
electron (hole) momentum at energy E hy k (q), and 
setting |AOi3| = /, |i?(5A| = I, we can write for the elec- 
tron (hole) wave function the following expressions for 
the phase gains along trajectories: 

AOB 

MA) ^ Ve(^)e*(^'+""\ 

^OB 

M^} ^ ^/.(A)e^(9'-™), 

. _ AQB [') 

, - _^ AQB 

V'h(B)e*(9'^™) MB). 

Due to the proximity effect, the quasiparticle wave 
functions in the normal region are composed of electron 
and hole components, which are mixed by Andreev re- 
fiections on NS boundaries. Thus, we have 

*AOB (x) = (x) + bj;t^ (x) , (8) 

"fBQAix) ^aij^{x)+hp;^{x). (9) 

Here, the arrows indicate the direction of the group 
velocity projection on the x-axis. The amplitudes of 
normal and Andreev reflection are denoted by TZ and 
A = — «|y^|e'''*'^, respectively. 

For the phase of the latter we have taken the result that 
is valid for a clean NS boundary and for a quasiparticle 
exactly on the Fermi surface (^ = E' — Ep = 0). The 
unitarity conditions, which encode the notion that, e.g., 
an electron propagating along BQA must undergo either 
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Andreev or normal reflection at A, can be written as: 

aV'e ( A)e*('='+"^") = TZa^Je (B) - i |^| e~'^« b^Jh ( A)e*(«'-""') , 
HhiB) = nbi:hiA)e'^'^^-'""^ ~ ^|y^|e*x«a^/'e(B). (10) 

In the next section, we shafl use Eq. Hl()|l to find the 
density of Josephson current-carrying levels, which will 



yield the Josephson currenliii. 

II. DISPERSION LAW 

In the above system of equations, we can (up to 
normalization and immaterial phase factors) assume 
V'e(^) = i^h{A) = iPeiB) = V^ft(B) ~ 1, which leaves 
us with a homogeneous system of linear equations for 
the quantities {a, a, 6, b}, with matrix 



M = 



n 


^i{kl-{-7rw) 







-i\A\e'^^ 







7^gi(<?'-'n'"') 


^i{kl -\-7rw) 


7^ 










~i\A\e'^'' 




-1 



(11) 



The dispersion law for the supercurrent-carrying quasi- 
particles follows from the solvability condition, det M = 
0. Henceforth, we explicitly take / = I, which is true in 
case of specular normal reflection and a rectangular nor- 
mal region. Introducing x = Xb ^ Xa, we have 

detM = 2e*('=+'?)' -i:*, (12) 
D = cos {uTT + l{k - q)) + \A\^ cos (x - n-n) 

-\n\^cos{k + q)l. (13) 

At this point we proceed by making an assumption 
about the energy-dependent amplitudes for Andreev and 
normal reflection. As we are interested in the possibility 
of normal reflection at the interface, we take these from 
the results of BTK. Thus, at low energies, we have& 



A 



n 



(14) 
(15) 



The momentum of an electron (hole) with energy ^ 
is fc(±^); thus we introduce the dimensionless quan- 
tity 95(0 = fc(^)/:/2 ^ ^2m{EF+i) ■ I and define 
■^±(0 = <^(0 ± <^(-C), giving 



-C(^, v) — cos((/)_ -f TTV) -f 
A2 



4Z2(A2 - e) 

cos 0+ — COS . 



cos(x — Tlil) 
(16) 



In order to calculate the density of states, we will 
use the approach due to Slutskini^ (applied to SNS sys- 
tems in Ref . IT^ to extract the density of current-carrying 
states iV(^) directly from the dispersion law. Thus, 
iV(^) = (5(5 — 5s), where 5s are given by the solutions 
of dispersion equation D{S^,u) = 0. A simple variable 



change yields 



95 



az?(5,^) 



95 



5(i?(5,z^)) 



21 . . 
sm( 



(17) 
(18) 



where in Eq. p8|) we have displayed the low-energy limit 
of 9-D(5, i')/95 and approximated 90_/95 ^ 21 /vp and 

a0+/95 - 0. 

The delta-function constraint in Eq. H17() implies that 



± arccos 



{2m 



A2 



V4Z2(A2 -52- 



(19) 

[cos(x - nQ) + 1] - 1 



with m being an integer. Let us briefly pause to discuss 
(qualitatively) the physical meaning of Eq. H19|l . The 
argument of the arccos in Eq. (|19(l has the interpreta- 
tion of a generalized phase difference across the junction. 
Since the current is given by differentiating the free en- 
ergy with respect to the phase differenceii, we see that 
the ± in Eq. p9|) correspond to levels that carry cur- 
rent in opposite directions. We shall refer to these two 
distinct contributions to N{£^) as -/V±(5)- The prefactor 



\dD/d^\^^\sm{ 



TTi>)\ assigns a particular weight 



to the (5- function contributions to N±; a direct applica- 
tion of the Poisson summation formula^"* yields (taking 
the low-energy limit) 



(20) 



T{x) = arccos [cos x + 1] - 1 ) , (21) 

where we have have rescaled the barrier height Z via 
Z ZA and we choose a branch of the arccos in Eq. H21|) 
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such that r smoothly goes to x — Trfi for Z 1/2. In 
the next section, we apply Eq. (|21|) to the calculation of 
the current due to a single trajectory, which (as we have 
briefly noted above) is related to the difference iV+(C) ~ 



it does not include, e.g., non- horizontal straight trajecto- 
ries, unless they are a part of a (closed) class A trajectory. 

IV. SUMMATION OVER CLASSES OF 
TRAJECTORIES 



III. JOSEPHSON CURRENT DUE TO A 
PARTICULAR TRAJECTORY 

Having determined the density of states A^(C) associ- 
ated with a particular trajectory, we now turn to the com- 
putation of the associated current. In the limit of zero 
normal reflectance the two terms A^±(C) in Eq. H21(l cor- 
respond to two separate groups of Andreev levels (AOB 
and BQA) which carry current in opposite directions and 
also move in opposite directions as x is changed^'"'. The 
simplicity of this expression is deceptive: the length / and 
phase gains f2 depend nontrivially on the shape of the 
trajectory. The associated contribution to the Josephson 
current is given in terms of N±(^) asi^ 



It, 



- I dfternh ( — 
W-oo V 2 , 

(22) 

where e is the electric charge; henceforth we shall set 
e = 1. In Eq. if^ . the subscript "tr." denotes that this is 
the contribution due to one particular trajectory, such as 
that pictured in Fig.^ At low energies, ip-iO ~ 2^1 /vp. 
Utilizing the integral formula 



tanh 



2m (1 - 5nfi) 
/3sinh(27r/n/ (/3?;f))' 



we find for the current 



) cos {-Kvn) VFx 



-J f3vF sinh (27rZn/ (/3uf)) 



(-!)"+! sin nr(x-7rl7). 

(23) 

Equation together with Eas. lf^ . (|^ . is the central 
result of this paper—. 

For closed-loop trajectories we must take into account 
that there always will be a related contribution from 
the same trajectory, but in the counterclockwise direc- 
tion (we completely neglect the dynamical effects of the 
magnetic field) . In the next section, we will consider the 
dominant classes of trajectory which contribute to the 
current. Before doing so, however, we pause to note that 
in the limit ^ 0, Z ^ 1/2, Eq. reduces to the 
Ishii's formula* for a wide clean SNS junction. 



(0) 



8vFx 



-(-l)"+isinnx. 



/3uf sinh (27r/n/ (Pvp)) 

(24) 

Despite this similarity, our expression Eq. (|23|1 gener- 
alizes Eq. (|24|) by incorporating normal reflection. In 
particular, our expression is not expected to yield a 2<I>o- 
periodic pattern seen in the theoretical results'*,-^ because 



In the present section we apply the results of the previ- 
ous section to calculate the contribution to the Josephson 
current due to various types of quasiparticle trajectories. 
We assume that each such contribution has the form of 
the general expression Eq. (|23|l but with values oiv, I and 
51 taken from the geometrical properties of the associated 
trajectory. 
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FIG. 2: Typical straight trajectory contributing to the 
Josephson current. The intercept of the trajectory along the 
segment CF is the coordinate y in Eq. I|25|l . 

We begin by discussing the leading-order contribution 
due to straight trajectories, which undergo no re- 
flections from the sides of the device. Due to the fact 
that normal reflection is a specular process, the trajec- 
tories must have incoming angle = 0, as depicted in 
Fig. 121 It is easily seen that for them, the phase gains 
[i.e. Eqs. 14I6|) ] due to the magnetic fleld are 1^ = and 
Q, = —2y(/)/W, where = ^/^o with $0 the flux quan- 
tum. The length I is given by L, and vpx — vp. Such 
trajectories are labeled by the coordinate y, given by their 
intercept on the line segment Ci^ in Fig. [3 To incorpo- 
rate all such trajectories, we integrate over y from — IF/2 
to W/2, dividing by an overall normalization factor W. 
Thus, we obtain for the current Is due to the straight 
trajectories 



Ltt 



dy 



00 ^ 



l)"+i 



^-^ sinh 1 



sinnT{x-2y(f>), (25) 



where for simplicity we have rescaled the dimensional 
coordinate y — > Wy to make it dimensionless and have 
defined w = W/L to be the dimensionless width of our 
sample. The parameter x = 2'kL/v-p[3 is a dimensionless 
measure of the temperature. 
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FIG. 3: A typical bouncing trajectory. The intercept of the 
trajectory on the hne segment CF is the coordinate y in 
Eq. 



Now consider the contribution due to class A (closed) 
bouncing trajectories, which bounce specularly once 
from the walls of the device. The requirement that the 
trajectory bounce specularly means, that there is only 
one for each intercept y. To compute this contribution, 
we must calculate the values of Z, and Q. These are 
found geometrically to be 



I = L^/l + w^ 
V = 2</>(y2 _ 1) 



0, 



(26) 
(27) 
(28) 



where we have again rescaled y — > yW . Summing over 
all such trajectories amounts to integrating over y (and 
normalizing by an overall factor of XjW as before) and 
multiplying by an overall factor of 2 to include counter- 
clockwise trajectories. Thus, we find for the current /b 
due to the bouncing trajectories: 



dy 



cos(7rz^n) 



sinh nx\f\ 



X (-!)"+! sin 7ir(x). 



(29) 



The generalization of Eq. H29|) to the case of multiple 
reflections from the sides of the device is straighforward. 



V. RESULTS AND DISCUSSION 

Having noted the various dominant contributions to 
the current, we now turn to the calculation of the critical 
Josephson current /c as a function of the applied flux. 
This quantity is defined in terms of the total Josephson 
current / (= /g + /b) as 



We begin by noting that the contribution due to the 
bouncing trajectories (i.e, Eq. 123) is explicitly depen- 
dent on the width-to-length ratio of the junction w and 
indeed vanishes for large w. Thus, we can effectively iso- 
late the effect of the straight trajectories by examining 
a very wide sample. In Fig. ^ we display numerical cal- 
culations of /c(0) for the case w = 50, Z = 1,2, 10, and 
x = (i.e., T = 0). (All curves have been normalized to 
their values at (/) = 0.) 




/c(0) EE max{/(x,0)}|o< 



<X<27r- 



(30) 



FIG. 4: The evolution of the normalized critical current (com- 
puted numerically) as a function of the external flux ^/"I'o for 
the case of a wide junction. The dot-dashed curve is Z = 1, 
the solid curve is Z = 2 and the dashed curve is Z = 10. 

The pattern in Fig. 01 is <i>o-periodic (i.e., Ic first van- 
ishes at (f> = ±1), which is expected for a wide junction^. 
For a barrier height on the order of the pair potential 
(Z = 1, dot-dashed curve), there is clearly a flattening 
of If. below (f> = ±1 arising from the energy-dependent 
amplitude for Andreev reflection. This flattening (but 
not the sharp peak for ~ 0) is reminiscent of the flat- 
tening effect seen in the JHAT experiment''. This ef- 
fect disappears for larger barrier heights, as exhibited by 
the Z = 10 curve (dashed line) in Fig. ^ This shows 
approximate reverting to the Fraunhofer shape is to be 
expected, as the contribution of higher-order processes 
to the Josephson current is suppressed by increased nor- 
mal reflection probability. At smaller values of w, the 
increased importance of the bouncing trajectories leads 
to somewhat different behavior of the Ic{(t>) curve as a 
function of Z. In Fig. |S1 we display /c(0) for the case of 
a narrow junction (i.e., w = I) for the cases Z = 1, 2, 10 
and 101 < 1. This curve is also for the case of T = a; = 0. 
The Z — 1 curve (dashed line) does exhibit flattening for 
small values of (f). In contrast to the wide junctions, in- 
creasing the barrier height leads to a sharpening of the Ic 
curve near (f> = 0, as seen in the Z = 10 curve (dash-dot 
line). 



VI. CONCLUDING REMARKS 

In this paper we have studied the effects of a finite tun- 
neling barriers and restricted geometries on the Joseph- 
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FIG. 5: The evolution of the normalized critical current (com- 
puted numerically) as a function of the external flux $/$o for 
the case of a narrow junction (w = 1), indicating a progres- 
sive sharpening of the Fraunhofer diffraction pattern as Z is 
increased. The dot-dashed curve is Z = 10, the solid curve is 
Z = 2 and the dashed curve is Z = 1. 

son effect in a clean junction. We derived closed formu- 
las for the density of states and Josephson current in the 
system, which generalize the known results to the case 
of BTK normal scattering on NS interfaces, which allow 
an intuitive insight in the mechanisms of non-Fraunhofer 
diffraction patterns in mesoscopic SNS junctions and can 
be used in complex geometries. 

As an application of the technique, we have calculated 



the critical Josephson current Ic in a rectangular SNS 
junction for various values of the junction width and the 
tunneling barrier height Z. As expected, the deviation of 
Ic{4>) from the case of a clean wide junction is most pro- 
nounced at small 4>, leading to a flattening effect, reminis- 
cent of the experimental results of JHAT— . Nevertheless 
there remain open questions. 

First, our calculations do not provide a quantitative 
match. More importantly, they do not reproduce the 
fine oscillations superimposed on the flattened first peak, 
seen in the experiment of Ref. 0. A possible source of 
this effect can be due to the closed trajectories with mul- 
tiple reflection points, not taken into account here. The 
fine structure in the density of states in SN billiard struc- 
tures (i.e., so-called Andreev billiards) arises from such 
trajectories^, in accordance with the Gutzwiller formula, 
which make this scenario at least plausible. Finally, we 
neglected the elastic scattering inside the normal region, 
which was significant in Ref. vl. Generally, it would keep 
I I, so that that calculations up to Eq. would hold, 
but the following formulas would be modified. 
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